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I. INTRODUCTION 

Recent developments in the application of variational 
methods [IHS] to large bases of hadron interpolating fields 
have made the extraction of the excited spectrum of 
hadronic states a realistic possibility (see e.g. 0HS])- 
Moreover, these new developments in lattice QCD pro- 
vide a solid foundation for studying multi-hadron scat- 
tering in lattice QCD [6]-|T2]. Because lattice QCD is 
formulated in Euclidean space, we do not have direct 
access to scattering amplitudes [12] • Fortunately, in a fi- 
nite volume, Luscher formula |14l presents a way to relate 
the discrete energy spectrum in lattice QCD to scatter- 
ing amplitudes. Since then, the framework derived by 
Luscher in [Tl] has been extended to moving frames [T51 - 
[19] , and the inelastic region by including coupled-channel 
effects [MIS]. 

Based on the approach developed in [25] for the scat- 
tering of two-particle system in finite volume, in this 
work, we discuss the three-particle system, consider- 
ing the finite-volume representation of the isobar model 
[2"oT, |2"7] in which interactions in the three-particle sys- 
tem are approximated by two-body scattering. We de- 
rive expressions which relate finite-volume energy shifts 
to isobar-model production amplitudes. 

The paper is organized as follows. In Section |ll] we 
discuss the three-particle system in finite volume. The 
summary and outlook are given in Section [III| 



II. THREE-PARTICLE SCATTERING ON A 
TORUS WITHIN AN ISOBAR APPROXIMATION 

In QCD, because of the low mass of pions, three- 
particle states become important even at relatively low 
energies. Methods for extracting three-particle ampli- 
tudes from lattice simulations are thus needed. In gen- 
eral, the scattering of a multiple-particle system is quite 
complicated, however, in some cases the quasi- two body 
approximation has proven to be quite successful [25] HH] • 
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A recent consideration of the three-particle system on a 
torus in the framework of Faddeev equations is presented 
in [3D]- 

In this work, we are interested in establishing the pres- 
ence of three-particle and two-particle resonances in a 
three-particle final-state system. In lattice QCD, the 
multi-hadron interpolating fields operators O are related 
to the total energies of multi-hadron system through the 
correlation function C(t) = (0|O(t)O t (0)|0) by relation 

C(t) = ^(0|O|n)(n|Ot|0) e - £ " t . (1) 

n 

Therefore, instead of building the relations between 
three-to-three particles scattering amplitudes and E n in 
finite volume in [31j . we intend to reconstruct a three- 
body resonance by its decay products, see figure [l] When 
two-body resonances occur, the decay or production of 
a three-particle final state system can be approximately 
described in terms of isobars recoiling from a bachelor 
particle [M [37J 112 As an example, the decay pro- 
cess of the strange axial meson resonance i^i(1400) to 
the Ktttt final state can be quite well approximated by 
if* (892) 7T in 5-wave, with the K* isobar decaying to Kit 
in a P-wave. 

As presented in Appendix [S] we consider the interme- 
diate resonance to be effectively an asymptotic state and 
parameterize its decay via a two-body isobar resonance 
recoiling from a spectator particle. 

In this section, we will generalize the methods devel- 
oped in [25] to apply to the non-relativistic three-particle 
system in the isobar approximation. For simplicity we 
consider three scalar particles with equal masses, m. As 
discussed in previous paragraph, we consider the pro- 
duction of a three-particle final state from the decay of 
a sharp resonance, \R), which we will treat as an asymp- 
totic state whose origin is in the short-distance qq dy- 
namics, rather than the final-state interactions. We pro- 
ceed assuming that the dynamics of transition from a 
resonance state in the far past to a three-particle state 
in the far future is described by a three-particle final- 
state interaction Hamiltonian H FS1 = H° + V FSl . Under 
the assumption of strong pair- wise final state interactions 
between two particles, strong enough to generate a two- 
body isobar resonance, the final state interaction Hamil- 
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FIG. 1: A schematic view of the isobar approximation to three-particle scattering, (a) The general 3 — > 3 amplitude is 
assumed to be dominated by two-particle 'isobar' resonances and to pass through a sharp three-particle resonance, R. (b) We 
assume the relevant singularities of the scattering amplitude are captured by considering the decay of the sharp intermediate 



resonance. 



tonian of the three-particle system takes the form 



H 



FSI 



fc=l,2,3 



p1 

2m 



fc=l,2,3 



where (pfc,Xfc) are the momentum and coordinate of the 
k th particle. 

The amplitude to produce the three-particle final state 
IP1P2P3) is denned by 



and the isobar. The solution of the Schrodinger equation 
H FS1 tfjjM = E ipjM is given by 

4 j M(nj, r*) - J d% d 3 v' k Goirij - r> ij} r fc - r' k ; E) 

x l> x 2i x 3)' 

with Go being the non-relativistic three-body Green's 
function, 



/(Pip 2 p 3 ^i?) = (PiP2P 3 |^ bl |* m (0)>, (2) G (v lJ ,v k ;E) = 



where the incoming wave state at t = is 



given by |* m (0)) = U(0, -00] 



exp 



R) and U(t,t ) = 
> . The wavefunction 



-i \t diV ffOt V FSI e-^ 0i 
j to 

of three final state particles in coordinate space is given 
by 



V>jm(xi,x 2 ,x 3 ) = (xix 2 x 3 |* in (0)). 



(3) 



Under the isobar model approximation, where V FS1 in- 
cludes only pair-wise interactions and where rearrange- 
ment among different isobar pairs is ignored, the three- 
particle production amplitude is approximated by 

/(P1P2P3 <- R) « J2 (^jk k \t^\R), (4) 

fc=l,2,3 

where is the two-body T-matrix operator for pair 
(ij), the relative momenta (qij , k^) are given by qjj = 
|(Pi - Pj) and k fc = |(P» + Pj - 2 Pfc) respectively. The 
S'-matrix describing the transition from \R) to IP1P2P3) 
is discussed further in Appendix |A"| 

In the center-of-mass frame, the wavefunction of the 
three-particle system has the following form [2"5] 

■0jm(xi,x 2 ,x 3 ) = ^2 ^JM( r «' r *)) 
fc=l,2,3 

where = Xj — Xj is the relative position between the 
two particles forming the (ij) isobar and = |(xj + 
x,) — Xfc is the relative position of the spectator particle 



d 3 q lj d 3 k k 



oiqyry+ikfc-r* 



(2tt) 3 (2tt) 3 



i? — 3m 



_ 3k| 

m 4m 



(5) 



Using the identity 



1 



E - 3m 



_ 3k| 

rn 4m 

1 



1 

2^ 



3k 2 n 2 

E-m- - ^ + zO S y - 2m - ^ + zO 

(6) 

We may represent the three-body Green's function by a 
product of two independent two-body Green's functions 



Go (r (j ,r fe ; E) = ^ / dEy G^ )+fc (r fe ; £) G +J (r 2J ; E l3 ) , 



where 



G 



(ij)+fc 




(r fe ;£) 



oikfc-I-fe 



d 3 k fc 



(2tt) 3 E-m-R 



_ 3k| 
8 J 4m 



is the free Green's function for propagation of the system, 
of total energy E, made of the isobar pair (ij) (with 
invariant mass Eij) plus the spectator particle, and where 



rf 3 q 2 



( 2 -) 3 ^ - 2m - < 



is the free Green's function describing the propagation of 
particles, inside the isobar. 
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The Green's functions are explicitly given by 



2tt 



E ~ m Am 2 e lqE ^ e lkE ^ rk 
dE *J ^7T^2 



2)ii 



3(4tt) 2 



r 



1.1 



(7) 



G {vij,r k ;E) -> 
with 

q Etj 
k Eij 

In Eq.|7), J^^dEij has been replaced by m dEij, 
becausexhe three-body Green's function has oscillatory 
behavior for Eij only inside the physical region. The 
three-particle wavefunction is given by 



m(Eij — 2m), 
-m(E — m — Eij). 



V>jm(xi,x 2 ,x 3 ) 
i \ - 



2tt ^ 

fe=l,2,3 



L k M L . 



2 m 



x frsy (lE zj n 3 )i s ^ Y s%j m Sxj (ry) 

x ht k {k Ei] r k ) i Lk Y LkMLk (r k ) 

X fjM;S tj M Stj ;L k M Lk (lEy > k E zj ), 

where the three-body production amplitudes are defined 
by 



/ 



JM:,S li ,Ms ij ;L k M Lk 

Am 2 



3 

xr^ jLt (fc Ej/ [)y^ fc (f' fc ) 

x ^i(r^)^M(xi,x^,x' 3 ). 



(8) 



The total angular momentum of the three-particle sys- 
tem is J, made up of (/Sy, Ms id ), the spin and its z-axis 
projection of the isobar pair (zj), and (L k ,M Lk ), the rel- 
ative orbital angular momentum and projection between 
fc th particle and isobar pair (ij). 



The partial wave production amplitude can be 
parametrized by 



fto) 

3 JM-S %i M s ,-L k M L , 



= fa M Sj . ; L fc M Lfc | JM) / Sij (q Ei . ) / Si . Lfc , 7 (fe Sy ) , 



(9) 



where 



, / x 47T iS 

/s„(<tey) = — e s 



sin 5 S 



is the partial wave scattering amplitude of two spinless 
particles inside the isobar pair (ij) and fs i:j L k j(k Eij ) is 
the partial wave production amplitude of the isobar pair 
(ij) and the spectator particle. 



The final expression for the three-particle production wavefunction, including the homogeneous term, reads 



^jm(xi,x 2 ,x 3 ) 



I 

2n 



k=l,2,3S ij Ms i 
L k M Ls 



{S l3 M Sij ; L k M Lk \JM) Y SijM . (hj) Y L „ M . (r k ) 



E-r 



dEij I 



(4tt) j Sij {q Eij Tij ) + i q ElJ (q Ei . ry ) f Sij (q Eij ) 



■'/.'., h t,{k Elj r k ) fs^L k j(k EtJ )] 



- 1 [( 4?r ) 3h k (k Etj r k ) + i k Eij „ Lk 

Boosting the three-particle system from the center-of-mass frame to a lab frame with the center-of-mass of isobar pair 
(ij) fixed at the origin, Xj + =0, the wavefunction of the three-particle system at the lab frame can be written 
as the product of a plane-wave, e lP -( x i+ x 2+x3)/3^ anc j a pj ece depending only on relative coordinates, ry = Xj — Xj 
and r k — —x k . As in the case of two-particle scattering |15) , requiring periodicity of the lab frame wavefunction with 
respect to ry and , we find that the boundary condition of CM frame wavefunction reads 

V^(ry + n^L, r k + n k L) = e^ n * L ^(Wk). 
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where k = 1,2,3 and Q is the Bloch wave vector for the three-particle system. The connection of the Bloch wave 
vector Q to the total momentum of the three-particle system is given by P = 37Q. Using the periodicity of the 
potential V^(r^ + n^L) — V i3 -(r- ■), the three-particle Lippmann-Schwinger equation on a torus can be written 



i>jM • (xi,x 2 ,x 3 )= E 



L 3 



fc=l,2,3 ' 

where the periodic three-body Green's function is given by 



G Q {r ij ,r k ;E)= Go(r tJ ~ n, tJ L,r k - n k L; E) e l ^ n - 



or asymptotically 



G Q (rij,r k ;E) -> 



2n 



dEa—^- — -2 2 T3 / , 



„ik k -r k 



3 L ^ <-^ L ^ Q ^-^ 



where P Q = {q G M 3 |q = + Q, for n G Z 3 }. Next, we use the expansion of the Green's function in Eq.(B2) in 
[2"5] and the definition of the scattering amplitude in Eq.Q and ([9|, to obtain 

^ L M Q) ( Xl ,x 2 ,x 3 )-> ^ E E E (^■A^;i fc M i jjM)y^. M ^(f iJ )r i , M , fc (f fc ) 

fe=l,2,3SijM Sjj . S'^M's.. 

L k M L k L'M'j 



E-m 
X / dEi 

1m 



'IK fSij ('//•:. 



S StJ m Sz] ,3'jM's.. n s <. {q Eij r w ) - M^,., iS , (g^ ) Js' zj r y ) 



S Lk M Lk , L > k M> Lk n L > k (k Eij r k ) - M ( ? k \ lLk L , kM , (k Eij )jL' h (k Eij r k ) 



Matching a general wavefunction of form J2jm c JM ^m( x 1i x 2, X3) with X)ja/ V'xjwr ( x i, x 2 , X3), and pro- 
jecting out the partial waves, neglecting rearrangement effects from crossed channels as appropriate in an isobar 
approximation, we get a set of equations which have general form m (iEy F(Eij, ry, r^) = 0. Since the variables 
( r iji r k) can be chosen arbitrarily, F(Eij,rij,r k ) = must be satisfied for each E^. From this we obtain three 
determinant conditions, the first is 



dot 



S Sij M 3iJ ,S> 3 M> Sij COt § Sij (q Eij ) - -MgyMsy ,S^M^ (fei ) 



= 0, 



(10) 



which is Luscher's formula for scattering between i th and j th particles inside the isobar pair (ij). It provides the 
constraint on the phase shifts 5s i:j as a function of invariant mass of the isobar pair (ij), Eij. The second condition, 



det 



^rfs t] L k j{k Eij ) 



+ i 



(SijM Sij ; L' k M' Lk \J'M') (SijM Sij ; L k M Lk \JM) M 



(Q) 

L k M L ,^LiM 



M Sij M H M' Lh 



k • k Li, 



(11) 



is a generalized Luscher's formula in moving frames for production of the spectator, k th particle, and the isobar (ij) 
with the specific spin Sij and mass Eij from initial state with spin J. Thus, Eq. (11) gives the constraint on the 
production amplitudes fs i:i L k j as function of the total energy, E, for each individual partial wave of isobar pair (ij), 
. 

The final condition, 



det 



SjM.J'M' S Slj ,S' zj h k ,L' k cot8 Sij (q Eij ) 



^rfs tj L k j(k Eij ) 



+ i 



J2 (SijM' Si/ ,L' k M' Lk \J'M') (S i jM Sii ,L k M Lh \JM)M { °l MSijiS ,^^ (k Eij ) 



= 0, 



(12) 
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leads to an additional constraint on both 5,g.. and fs i3 L k j for the scattering between the k th particle and all the 
allowed partial waves of the isobar (ij). 

I 



In the case that only a single partial wave Sij of the 



isobar pair is dominant, Eq. ( 12 ) becomes redundant and 



the three conditions reduce to two, 



cot 5 Sij {q Eij ) = M 



(o) / \ 

and Eq. ( |TT| ). Because of a finite- volume with cubic 
boundaries, the continuous rotation symmetry is reduced 
to the little group of allowed cubic rotations that leave 
the centre-of-mass momentum invariant, so that, Eq. 
(10 1, (11) and (12) have to be subduced according to ir- 



reducible representations of the appropriate little groups 

M- _ _ _ 

From Eq. (10), (11) and (12), we find that even in 
relatively simple cases, e.g. a single isobar dominating a 
single relevant partial wave, extracting the phase-shifts 
and determining the invariant mass from the mea- 
sured total three-particle energy E is a rather difficult 
task. For additional information, we can first perform 
computations of two-particle correlators with the quan- 
tum numbers of the isobar channel to obtain information 
on Eij. However, in the three-body calculation, there 
may be multiple E^ for each individual E allowed by 
kinematics (2m < E^ < E — to), thus, finding the cor- 
respondence between Eij and E may requires some as- 
sumptions and model input. 

In recent works [U \5\ , determining the spin of the ex- 
cited states by considering the overlap of carefully con- 
structed operators with the particular state, (n|O|0), has 
been proven to be successful in certain cases. We may 
use a similar idea to identify the E^ and E relation 
in a three-particle system, by considering the overlap 
of operators with the state having the particular quan- 
tum numbers (Sij , , J) and invariant mass of the isobar 
pair E^. For instance, working in center-of-mass frame, 
the three-particle operator may be constructed in such 
way that an isobar pair (ij) has definite relative momen- 
tum, ||pi — pj| = ^|n g |,n g G Z 3 , and definite total 
momentum, |p, + pj\ — ^|nfe|,rifc e Z 3 . This opera- 
tor will strongly overlap with the state having invariant 
mass of the isobar pair (ij), Eij ~ 2m 2 + ^(^-n q ) 



amplitudes can be factorized as the product of two indi- 
vidual scattering amplitudes, one describes the scattering 
of two particles inside isobar pair, and another describes 
the scattering between isobar pair and the spectator par- 
ticle. Three determinantal conditions, Eq. (10), (11) 



and (12), were obtained for three-particle scattering in 
a finite volume. One condition, Eq. (10), relates the 



scattering phase shifts of two particles inside an isobar 
pair to the invariant mass of the isobar pair. The other 



two conditions, Eqs. (11) and (12), relate the scattering 



phase-shifts between an isobar pair and the spectator to 
the total energy of the three-particle system. A proposal 
for extracting the phase-shifts of a three-particle system 
from lattice QCD simulations is presented that makes use 
of carefully constructed two-particle operators within the 
overall three-particle operator construction. 
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Appendix A: Production amplitudes and Faddeev 
equations 

Our task here is to consider the transition from a res- 
onance, treated as though it were an asymptotic state, 
in the far past, to a system of three particles in the far 
future. In this case the S'-matrix [35] for the transition 
reads 



(p lP2 p 3 |C/(oo,0)C/(0,-oo)| J R> = (*° ut (0)|v]/ in (0)>, 



(Al) 



and the total energy of the three-particle system E ~ where U(t,t ) — Tlexp — i f* df 'e lH<>t V FSI e lH " t ' } 



to + E^ + -^(^-iik) with small shifts caused by the in- 
teraction. Additionally, in the case of isobar dominance, 
fermion bilinear operators subduced from spin Sij are 
likely to have good overlap. 



III. SUMMARY 

For the case of three non-relativisitic particles undergo- 
ing scattering, we considered in finite-volume the isobar 
model approximation, where quasi-two-body scattering 
is assumed to be dominant. Under the isobar model ap- 
proximation, the three-particle partial wave scattering 



and the incoming and outgoing wave states at t = are 
defined by 

|* in (0)) = U(0,-oo)\R), 
|M/° ut (0)) =£/(0,oo)|p lP2 p3). 

Standard manipulations lead to the production ampli- 
tude 

/(P1P2P3 <- R) = ( P iP 2 p 3 |U FSI |vl/ in (0)). (A2) 
The incoming wave state satisfies equation 



l* in (o)> 



1 



1 

E - H FSl + iO 



V 



FSI 



\R). (A3) 
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So that, the production amplitude is finally given by 

/(P1P2P3 <-R) = (P1P2P3I T FSI \R), (A4) 
where the T-matrix operator is defined by 



wise interactions, and thus the T-matrix for the three- 
particle system can be decomposed as T FSI = J2 k T^ k , 
where TW* satisfies Faddeev equations [25] . 



T FSI = V FS1 + V FSI 



1 



E - H FSl + te 



-V 



FSI 



1 



E - H° + ie 



(A5) 



Within the isobar model we assume dominance of pair- where t^' is the two-body T-matrix for pair (ij). 



It follows that the three-particle production amplitude can be decomposed 
/(P1P2P3 <— R) = J2 k f^' k (<lij kfc 4- R), where f^ k (<\ij k^ 4- R) satisfies an integral equation [3B] 



/fe) fc (q l3 k fe 4-i?) = 



+ 



<qyk fc |*«)|i*) 

d% k d% ( % k fc |^)|q^) fW*( q r k k<^R) 
(2^)3(2^)3 E-H^^ + ie 

d\' ki d% (q^kfc^|q^) /^(q'fcjkj 
(2^)3(27r)3 E-ff°(q' fci k'.) + ie 



(A6) 



as 



(A7) 



The first term on the right hand side in Eq. \K7\ is the 
Born-term for production of the isobar pair (ij) plus the 
k th particle from the resonance state \R). The remain- 
ing terms generate the rescattering effect from different 
isobar pairs - the isobar model approximation is to keep 
only the Born term in Eq. (A7), so that the production 



amplitude is approximated by 

/(P1P2P3 <- R) « Yl (<Hjkk\t (ij) \R). 

fc=l,2,3 



(A8) 



C. Michael, Nucl. Phys. B 259, 58 (1985). 

M. Luscher and U. Wolff, Nucl. Phys. B 339, 222 (1990). [16 

B. Blossier, M. Delia Morte, G. von Hippel, T. Mendes 

and R. Sommer, JHEP 0904, 094 (2009) [17 
J. J. Dudek et al. (Hadron Spectrum Collaboration), 
Phys. Rev. D 82, 034508 (2010). [18 
R. G. Edwards, J. J. Dudek, D. G. Richards, and S. J. 
Wallace, Phys. Rev. D 84, 074508 (2011). [19 
J. J. Dudek et al. (Hadron Spectrum Collaboration), 
Phys. Rev. D 86, 034031 (2012). [20 
X. Feng, K. Jansen, and D. B. Renner, Phys. Rev. D 83, [21 
094505 (2011). 

J. J. Dudek et al. (Hadron Spectrum Collaboration), [22 
Phys. Rev. D 83, 071504 (2011). 

S. R. Beane et al. (NPLQCD Collaboration), Phys. Rev. [23 
D 85, 034505 (2012). [24 

C. B. Lang, D. Mohler, S. Prelovsek and M. Vidmar, 
Phys. Rev. D 84, 054503 (2011). [25 
S. Aoki et al. [CS Collaboration], Phys. Rev. D 84, 
094505 (2011) [26 
J. J. Dudek, R. G. Edwards and C. E. Thomas, 
|arXiv:1212.0830l [hep-ph] . [27 
L. Maiani and M. Testa, Phys. Lett. B 245, 585 (1990). [28 
M. Luscher, Nucl. Phys. B 354, 531 (1991). 

K. Rummukainen, S. Gottlieb, Nucl. Phys. B 450, 397 



(1995). 

C.-J.D. Lin, G. Martinelli, C. T. Sachrajda and M. Testa, 
Nucl. Phys. B 619, 467 (2001). 

N. H. Christ, C. Kim and T.Yamazaki, Phys. Rev. D 72, 
114506 (2005). 

V. Bernard, Ulf-G. Meifiner and A. Rusetsky, Nucl. Phys. 
B 788, 1 (2008). 

V. Bernard, M. Lage, Ulf-G. Meifiner and A. Rusetsky, 
JHEP 0808, 024 (2008). 

S. He, X. Feng, C. Liu, JHEP 0507, 011 (2005). 

M. Doring, Ulf-G. Meifiner,E. Oset and A. Rusetsky, Eur. 

Phys. J. A 47, 139 (2011) 

S. Aoki et al. [HAL QCD Collaboration], Proc. Japan 
Acad. B 87, 509 (2011) 

R. A. Briceno and Z. Davoudi , [arXiv: 1 204. 1 1 10| [hep-lat] . 
M. T. Hansen and S. R. Sharpe, Phys. Rev. D 86, 016007 
(2012) 

P. Guo, J. Dudek, R. Edwards, A. P. Szczepaniak, 
|arXiv:1211.0929l [hep-lat]] . 

Y. Goradia and T. A. Lasinski, Phys. Rev. D 15, 220 
(1977). 

G. Ascoli and H. W. Wyld, Phys. Rev. D 12, 43 (1975). 
L. D. Faddeev, Mathematical Aspects of the Three-Body 
Problem in the Quantum Scattering Theory (Israel Pro- 
gram for Scientific Translation, Jerusalem, Israel, 1965). 



7 



[29] L. D. Faddeev, Zh. Eksp. Teor. Fiz. 39, 1459 (1960) 

[Sov. Phys.-JETP 12, 1014(1961)]. 
[30] K. Polejaeva and A. Rusetsky, Eur. Phys. J. A48, 67 

(2012). 

[31] R. A. Briceno and Z. Davoudi, [arXiv:1212.3398| [hep- 
lat]]. 

[32] R. C. Hwa, Phys. Rev. 130, 2580 (1963). 

[33] R. Aaron and R. D. Amado, Phys. Rev. Lett. 31, 1157 



(1973). 

[34] C. E. Thomas, R. G. Edwards and J. J. Dudek, Phys. 

Rev. D 85, 014507 (2012). 
[35] M. Gell-Mann and M. L. Goldberger, Phys. Rev. 91, 398 

(1953). 

[36] A. Ahmadzadeh and J. A. Tjon, Phys. Rev. 139, B1085 
(1965). 



